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Abstract 

The static solutions of the axially symmetric vacuum Einstein 
equations with a finite number of Relativistic Multipole Moments 
(RMM) are described by means of a function that can be written 
in the same analytic form as the Newtonian gravitational multipole 
potential. A family of so-called MSA (Multipole-Symmetry Adapted) 
coordinates are introduced and calculated at any multipole order to 
perform the transformation of the Weyl solutions. 

In analogy with a previous result pO] obtained in Newtonian grav- 
ity, the existence of a symmetry of a certain system of differential 
equations leading to the determination of that kind of multipole solu- 
tions in General Relativity is explored. The relationship between the 
existence of this kind of coordinate and the symmetries mentioned is 
proved for some cases, and the characterization of the MSA system of 
coordinates by means of this relationship is discussed. 
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1 Introduction 



As is known, the description of Newtonian Gravity (NG) in the vacuum in- 
volves solutions of the Laplace equation whose general well-behaved solution 
is a series with arbitrary constants that can be identified with the Multipole 
Moments (MM) of the source, and these quantities allow us to characterize 
the specific solutions given by the succession of partial sums of the series. 

In contrast, the static solutions of the axially symmetric Eisntein vacuum 
equations describing the gravitational field of a bounded isolated mass dis- 
tribution in General Relativity (GR) can be described by means of only one 
metric function, / = qqq, which satisfies the Ernst equation [T]. 

We are interested in the following questions: might it be possible to obtain 
a description of these solutions in GR by means of a function, namely u, with 
the same behaviour as the classical potential? Could we write this function 
u analytically equal to the Newtonian gravitational series but in terms of 
Relativistic Multipole Moments (RMM)? 

We are concerned with these questions for several reasons, in particular 
because such a description of the relativistic gravitational solution would re- 
cover the benefits of the classical interpretation of the gravitational potential 
(see pOl for details). Moreover, the Weyl family of solutions depends on 
arbitrary constants, a„, in principle without any physical criteria to choose 
one or another solution from them, whereas the function u would allow us to 
deal, in a very simple form, with the Relativistic Multipole Solutions. This 
has been the aim of some authors and their works devoted to obtaining so- 
lutions of the Einstein vacuum equations with a finite number of prescribed 
RMM. 

In this work we seek an answer to these questions by introducing a family 
of coordinate systems referred to as MSA (Multipole-Symmetry Adapted). 
The possibility of extrapolating the symmetries obtained in NG \lOl to GR, as 
well as characterizing the solutions with a finite number of RMM by means 
of group-invariant solutions, are the relevant features of these coordinate 
systems and the reason for their proposed name. 

In a work published recently [lOj, the existence of some kinds of sym- 
metries in NG has been proved, which makes it possible to extract from 
all solutions of the axially symmetric Laplace equation those with the pre- 
scribed Newtonian Multipole Moments. A family of vector fields that are 
the infinitesimal generators of certain one-parameter groups of transforma- 
tions can be constructed. These vector fields represent symmetries of certain 
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systems of differential equations whose group-invariant solutions turn out 
to be the family of axisymmetric potentials related to specific gravitational 
multipoles. 

By introducing these coordinates, the function u, which describes the 
static and axially symmetric vacuum solution with a finite number of RMM, 
should satisfy the same system of differential equations as the classical po- 
tential in NG, and the symmetries of these equations [ID] thus allow us to 
describe and determine the Multipole Solutions in GR analogously to the 
Newtonian case. 

Since the function u, which is transformed from ^foo, must fulfil the corre- 
sponding Ernst equation written in MSA coordinates, the following question 
arises: is it possible to obtain conditions on the change of coordinates, to 
choose the suitable gauge, from the extension of the symmetries to the corre- 
sponding Ernst equation? In other words, can the symmetry groups obtained 
for the Laplace equation and the supplementary equation [10] be extrapo- 
lated to the Ernst equation in that system of coordinates? And if so, could 
we establish theorems relating the existence of the symmetry of a system of 
differential equations to that system of coordinates? We shall see, at least 
for the Monopole case, that this relationship can be used to determine the 
MSA radial coordinate explicitly. 

We shall try to answer these questions along the work in the following 
way: 

In section 2, the MSA systems of coordinates are defined from the context 
of the multipole expansion of gravitation introduced by Thorne [2J, and the 
procedure to calculate these coordinates is shown for each set of multipole 
structure of the desired solution. The procedure first consists of performing 
a coordinate transformation from Weyl coordinates, preserving the Killing 
vectors and the asymptotically Cartesian behaviour. Second, we introduce a 
function u by redefining the goo metric component and we force this function 
M to be a solution of the corresponding Ernst equation written in the new 
system of coordinates. The results are addressed in Appendix B. Some com- 
ments about the behaviour and interpretation of the coordinates obtained 
complete this section. 

In section 3 we attempt to provide these coordinates with a meaningful 
interpretation by means of the existence of symmetries of certain differential 
equations. We recall that in these coordinates the static and axially sym- 
metric vacuum solutions with a finite number of RMM can be described as 
group-invariant solutions of the same system of differential equations that ad- 
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mits the symmetries obtained for the Newtonian case [TD]. Furthermore, we 
prove two theorems that extend those symmetries to the corresponding Ernst 
equation written in the MSA system of coordinates for the Monopole and the 
Monopole-Dipole cases. And these theorems provide a relationship between 
the existence of these symmetries and the Newtonian-type solutions in GR, 
at least for these two cases. In subsection 3.2 a possible characterization of 
these systems of coordinates for each multipolar solution is explored, without 
taking into account our knowledge of the corresponding set of constants an 
for each case. In this sense, for the Monopole case the corresponding MSA 
coordinates can be obtained as the unique solution of the Ernst equation 
(and the suitable constraints) with some boundary condition. Nevertheless, 
for any other case this procedure fails to provide the uniqueness of the MSA 
system, since two coordinates must be solved rather than the radial coordi- 
nate alone, as is the case only for Spherical symmetry. Finally, in Appendix 
A the expressions of the RMM in terms of the constants and the inverse 
relation are shown for a solution with a set of arbitrary RMM up to order 
10. 



2 The MSA system of coordinates 
2.1 Definition 

In 1980 Thorne introduced a system of coordinates called ACMC {Asymptot- 
ically Cartesian and Mass Centered) in the context of Multipole expansions 
of gravitational radiation [2]. His work presents a definition of Relativistic 
Multipole Moments (RMM) and shows us how to deduce the RMM of a 
source from the form of its stationary and asymptotically flat vacuum metric 
in an ACMC coordinate system. 

If the components of the metric are written in the coordinates {t, r, 6, 0} 
we can read off the RMM from the resulting expressions, other terms, R^ll~^^ (y) 
called Thorne rests, appearing at the same time that are functions depending 
on the angular variable y = cos 6 in at least one degree lower than those as- 
sociated with the RMM. For the case of axial symmetry, the qqq component 
of any static metric written in that kind of coordinates resembles: 



5'oo = -1 + ^ 



c2 



oo 



oo 

—,MnPniy) + j:7^Ao"\y) 



/y I L~\~ i- ^^^^^ nyii 

.ra=0 n=l 
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M„ being the RMM of order n, and P„(t^) the Legendre polynomiaL 

There is gauge freedom in the choice of ACMC coordinates preserving the 
invariance of the first series in ([1]) and addressing the differences in the metric 
expansion through the Thorne rests. Among the broad class of coordinate 
systems of this type, we wish that system to lead to an expansion of the 
metric in such a way that all the Rl^~^\y) Thorne rests will vanish. We 
propose that this system of coordinates should be referred to as ACMC-TRF 
{ACMC-Thorne Rest-Free) in a first step, and then become the so-called 
MSA {Multipole- Symmetry Adapted) system of coordinates for reasons we 
shall see in the next section. 

Thorne showed that de Bonder coordinates are ACMC-to order N for all 
N, and in [2] he discussed that a de Bonder transformation of the coordinates 
is not necessary. Hence, we do not impose the harmonic condition on our 
system of coordinates, and we perform a coordinate transformation of the 
following form 

r = x" + x"(xO ; x" = cte + 0(l//?) , (2) 

{x"} being the Weyl spherical coordinates {t,R,uj = cos9,{p}, and {x"} = 
{t,r,y = cos6,Lp} being the new system of coordinates. This keeps the 
Killing vectors unchanged and maintains the asymptotically flat form of the 
metric: ga/3 = rjafs + 0{1/R). Now, the x° functions must be adjusted so 
that the metric in the new system of coordinates will satisfy the following 
condition: (c = 1) 

^ 1 

goo = -l + 2u = -1 + 2^2 -^iMnPniy) , (3) 

n.=0 

where N stands for the number of RMM we wish to consider for our metric. 
Henceforth, the g^o component of the metric in this system of coordinates 
acquires a form that is analytically related to the Newtonian gravitational 
potential of a classical Multipole Solution, and ([3]) represents the solution, 
with a finite number (N+l) of RMM, of the static and axially symmetric 
Einstein vacuum equations. 
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2.2 Calculation of MSA coordinates for any Multipole 
Solution 



The hne element of a static vacuum metric is as follows: 

{i?, a; = cos9,(p} being the Weyl spherical coordinates, and \E', 7 are metric 
functions satisfying the following equation^ 



IR = ^ [R'^^l - (1 - uj')^l - 2i?CU^H^^] 

7^ = u;[R'^l-{l-uj^)^l]+2R{l-uj')^R^^ , (5) 

where the sub-indices denote partial derivation with respect to them. 

The general solution for an isolated source with axial symmetry is given 
by the following asymptotically flat series (the family of Weyl solutions): 

* = Eil^w. (6) 



n=0 



where coefficients arbitrary constants; any set of those coefficients 

univocally determines the solution. From the calculatioij^ of the RMM of 
this metric (jl]) with the function \I' one can obtain an expression for 
the coefficients a„ in terms of the RMM and hence it is possible to choose 
of these coefficients by neglecting the undesirable RMM. This procedure af- 
fords a Multipole Solution having a finite number of Multipole Moments [7] . 
Some authors have devoted some time to seeking such those solutions (Pure 
Multipole Solutions in GR) [7|, [S], [S]. In [7], the M-Q Solution is obtained, 



^Let us note that the integrabihty condition of the equation for 7 is simply the equation 
for ^E*, and therefore the solution of the Laplace equation univocally identifies the space- 
time. 

^We perform the FHP method [3], which allows us to obtain the RMM in terms of 
coefRcients m„ involved in the expansion series on the symmetry axis of a conformal 
Ernst potential. Since that conformal potential is related to the metric function 5* ([5]), 
the coefficients m„ can be expressed in terms of the set of coefficients and hence the 

final result provides a relation A/„ = M„(a„) (|54p : the triangular structure of this relation 
[2], [7j allows us to calculate the inverse relation, o„ = a„(M„) (155)1 . 
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as well as the Quadrupole Solution itself, and more recently in [TT], [12] a 

method has been proposed for obtaining the general terms of the series (the 

coefficients a„) that define the Pure 2^-pole Solutions. In fact, the general 

term of the series corresponding to the gravitational Dipole and the solutions 

with Monopole plus any other 2^-pole moment are written specifically. 

In Appendix A, the RMM of the general static solution (jH]) are written 

in terms of the coefficients well as the inverse relation for a solution 

with arbitrary RMM up to order 10. 
1 

Let us set M = -(1 + goo); this function u, corresponding to the solution 
([6]) written in Weyl coordinates (uw), resembles the following expression 



1 

Uw 



2 



n=0 



(7) 



and would provide the gravitational solution with a finite number {N + 1) of 
RMM by substituting the coefficients a„ from expression fl55|) . 

In terms of the Ernst potential [1], which is a real function for the static 
case, the goo metric component equals this potential, and hence the function 
u satisfies the following equation derived from the Ernst equation [T] : 

{2u-l)[R\RR + 2RuR + il-uj^)u^^-2uju^] = 2[R\l + {l-uj^)ul] . (8) 

It is straightforward to calculate the transformation of equation ([8]) by 
means of an arbitrary change of coordinates from the Weyl system {{R,uj}) 
to another one {{r,y}), leading to the following expressions: 

2 

AUrr + BUyy + 2CUry + DUr + EUy — ^ [^Mr + BUy + 2CUrUy] = 0, (9) 

where A, B, C, D and E are functions of the coordinates {r, ?/} defined as 
follows 

A{r,y) = LBi{r,r) B{r,y) = LBi{y,y) C{r,y) = LBi{r,y) 

D{r,y) = LB2{r) E{r,y) = LB^iy), (10) 

LBi{) and LB2O being the Laplace-Beltrami operators with respect to 3- 
dimensional Euclidean metric (with axial symmetry) written in Weyl spher- 
ical coordinates; i.e., 

LB,i ) ^ r/^^V.()V,() = R^dnOdai) + (1 - i^^)d^{)d^{) 
LB2O = v''^.,{) = R'dlj, + 2Rdn + il-uj')dl^-2cod^. (11) 
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Let us perform the above-mentioned coordinate transformation ([2]) by 
assuming the following asymptotically Cartesian behaviour of the new coor- 
dinates: 



r = R 



QO _ 

n=l 
oo ^ 

+ E^"H^- (12) 



oo 

y = w ■ 



n=l 

We shall now impose the following two conditions: 
(£"0) uw = umsa 

(EI) [AUrr + BUyy + 2CUry + DUr + EUy]^^^^^^ = 

^^[Aul + Bui + '^C^rUyl=u,^^^ ■ (13) 

The equation (EO) means that we force the solution with a finite number 
(A^ + 1) of RMM written in Weyl coordinates (uw) to be functionally equal 
to the Newtonian gravitational potential, with that number of Multipole 
Moments, written in the MSA system of coordinates (umsa) required to exist. 
This condition is equivalent to the coordinate transformation of the metric 
component goQ. The equation (EI) should be understood in the following way: 
the function umsa must be a solution of the differential equation ([9]) obtained 
from the Ernst equation for -u ([8]) by means of the gauge transformation. 

These two conditions univocally determine the functions fn{^) and gn{^^) 
(fT2|) . up to any order, in the following way. First, equation (EO) provides a 
relation between each fn{^^) and Qkij^), for k from 1 to n — 1, by developing a 
power series expansion on the inverse of the radial coordinate R. Second, we 
substitute umsa in the differential equation (EI) and the resulting expression 
can be expanded in power series of 1/R by using the gauge (IT^ . Since we 
have considered the previously obtained relations between functions /n(<^) 
and Qni^^) in this expansion, the condition (EI) finally leads to a complete 
determination of the coordinates {r,y}. 

We have computed the calculation of the gauge ( fT2l) for a solution with 
the first three RMM (Mo,Mi and M2, monopole, dipole and quadrupole 
moments respectively) up to order 0(l/i?^°); the results obtained are shown 
in Appendix B. 
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Let us offer some comments about the good behaviour of these coordinates 
{r,y}. According to the meaning of the RMM, these coordinates reveal the 
loss of relevance of high-order multipoles in the description of the solutions at 
large distances from the source. Moreover, as can be seen in the expressions 
of the functions Qni^^) fl58|) . all these functions vanish for u = ±1. This 
means that the coordinate y preserves the axial symmetry, since y = uj along 
the axis. They are not harmonic coordinates because that condition (the 
de Bonder gauge) is not fulfilled by the associated Cartesian coordinate 
In addition, if we consider the MSA coordinates for the case of spherical 
symmetry, which can be done by neglecting all rmM^I in (!56|) greater than 
the monopole (Mq), then we observe that all the functions /„(co' = ±1) = for 
n > 2, (i.e., the coordinate r along the axis orthogonal to the equatorial plane 
shows a good behaviour; in other words, it equals the Cartesian coordinate z, 
up to a displacement along the axis). In fact, the case of spherical symmetry 
deserves a more detailed analysis: the expressions (l57j) for this case can be 
written as follows: 

hniu;) = M^^Ci;y^\uj) , /2„+i(a;) = ,n > 1, h{uj) = M 

g2n{oo) = -M2"C^f(a;), w(^) = 0,n>l, g,{uj) = 0, (U) 

Cn"^^^"* being Gegenbauer orthogonal polynomials, and henceforth M = Mq. 
Since the generator function of these polynomials is known: 

oo 
n=0 

the following relations hold 

$:cL-^/^)(.)A- = l(r, + r_) 

n=0 

oo ^ 

-E^2n"l'^M^"^' = ^(-^---)' (16) 
n=0 

^Let {x^} = {x,y,z} be Cartesian coordinates associated with the spherical ones 
{f,9,(p} as z — fcos9, x + iy = fsin^e"^; these coordinates are said to be harmonic 
if □{i:*} = 0, where □ denotes the D'Alambert operator with respect to the metric con- 
sidered. This condition leads to the following equation for the coordinate z: LB2{z) — 0. 
It is easy to see that LB2{z) ~ cos 9LB2{f) + fLB2{cos9) + 2LBi{f, cos 9) and the coor- 
dinates {f = r, cos 9 = y} given in (fT2)) (|57ll58p do not satisfy this harmonic condition. 

''As can be seen, expressions ([56l) are recovered from ((57|) by taking all RMM, except 
for Mq, equal to zero. 
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and therefore, by taking A = M/R, the coordinates {r,y} from (fT^ for this 
case are given by 

R 

r = M+ — (r+ + r_) = M(x + 1) 
R 

y = ^^^+^^-'> = yp ^ (1'') 

where {x,yp} are the prolate spheroidal coordinates [18], [H]. 

This radial coordinate r is easily recognized since it is merely the so- 
called standard radial coordinate of Schwarzs child and, as is known, the 5^00 
component of the Schwarzschild metric written in this coordinates is goo = 
— 1 + 2M/r, and hence the prescribed form of the metric component ([3]) is 
recovered by this coordinate system and the relativistic Monopole Solution 
is described by a function u equal to the spherical Newtonian potential M/r. 

In the following section we shall see the reasons why we refer to this 
system of coordinates as a Multipole-SjTiimetry Adapted one. 



3 Interpretation and characterization of MSA 
coordinates 

The function u written in a MSA system of coordinates {r, y} should be a 
solution of equation (that is, the meaning of the condition (EI)), but 
at the same time it is also a solution of the following system of differential 
equations 

= r'^Urr + 2rUr + (1 — y'^)uyy — 2yUy = Au 

= d^+\, ' (18) 

for any value of N, whenever the function u represents the Multipole Solution 
with a finite number (A^ + 1) of RMM ([3]): 

N 

^MSA = Y.^My)- (19) 

n=0 

Therefore, the symmetries of the system of differential equations f|T8|) . 
obtained in NG [10], that allow one to extract from all solutions of the axially 
symmetric Laplace equation those with the prescribed Newtonian Multipole 
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Moments work identically in the case of using MSA coordinates, but now 
the quantities M„ are the RMM (1191) . In this sense, the same family of 
vector fields that are the infinitesimal generators of certain one-parameter 
groups of transformations can be constructed. These vector fields represent 
symmetries of the system of differential equations (fTSl) whose group-invariant 
solutions turn out to be the family of axisymmetric potentials related to 
specific gravitational multipoles ffT^ . 

The question we want to answer now is whether those groups of sym- 
metry exist for the system of differential equations joined by (fT8l) and ([9]); 
if so, we could generalize the above-mentioned results to GR, establishing a 
relationship between the existence of a certain symmetry and the uniqueness 
of the solutions of the Einstein equations with a prescribed multipole struc- 
ture. Let us remark that equation iQ is the corresponding Ernst equation 
for the function u, and hence we could say more appropriately that Einstein 
equations admit the symmetry. 

Moreover, we wish to know whether the action of these symmetries on 
the equation might provide conditions on the unknown functions /n(i^) 
and Qni^) ffT^ to characterize the MSA coordinates. 



3.1 Multipole symmetries in GR 

Let 

d d 

^ = r—-u— 20 

or ou 

be a vector field on an open subset M C X x f/, where X = is the 
space representing the independent variables, coordinates {x} = (r, y) being 
MSA coordinates such that (fT9l) is fulfilled for N = 0, and [/ = M, with the 
coordinate u that represents the dependent variable. 

We can state the following theorem: 
Theorem 1 

The system of equations Aj,(x,m*^"^) = given by 

Ai(x,m(")) = Au = 

A2(x, = = (21) 

A3(x,m(")) 

where Ai is the Laplace equation (with axial symmetry) ([IS]) , A2 the so-called 
supplementary equation [10], and A3 equation ([9]), admits a symmetry group 
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whose infinitesimal generator is v. 



Proof: 

Tlie prolongation of v acting on the supplementary equation is pr*^^W [A2] = 
—Uy, and the second prolongation of this vector acting on Ai is pr^'^^w [Ai] = 
— Am. Therefore, both prolongations are zero whenever the system of equa- 
tions A^(x,m(")) = is fulfilled. 

With respect to the third equation of the system, it is straightforward to 
see that 

pr^'^^V [A3] = DUr + EUy + r {ArUrr + Bj-Uyy + 2CrUry + D^Ur + Ej-Uy) + 

^2u~~l + ^'^yy + '^Cury + Dur + Euy) + 



2r 



2u - 



- [Aj-ul + BrU^y + 2CrUrUy) — 4 Ag 



(22) 



Since the prolongations of the vector only need to vanish on solutions of 
the system of equations (l2Til [13], we make use of equations A2 and A3 to 
obtain the following relation between the derivatives of u 



D 



Ur 



-U^ -Ur 



2m-1 " A 
and we substitute this and A2 into (!22|) to obtain 

D 



(23) 



pr'-^-'v [Ag 



rDr - —{rAr - A) 



+ ul 



2A 



{2u 



(24) 



The derivatives Ur and u^ are related by equations A3 and Ai (with Uy = 0) 
as follows 



Ur 



2A\ 



2A 



u^ 



r J ^ 2u — 1 
and therefore the second prolongation of v acting on A3 is 



(25) 



pr 



(2)v [As 



Ur 



(26) 



Equation fl23|) is the Ernst equation for the function u with the constraint 
A2. It should therefore be fulfilled by the Monopole Solution, represented by 
u = M/r, since we have required the system {r,y} to be MSA coordinates 
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for the Monopole case (spherical symmetry), and hence we can state the 
following relation between the coefficients A and D 

D 2 r-M 

— = • 27 

A rr-2M ^ ^ 

Finally, by using equation (!27l) and its derivative with respect to r, i.e.. 



in (126|) . we have that: 



,0) .AT 2AM 

pr^ '\ [A3J = :r^Ur 



1 



r{2u-l) r-2M 



(29) 



r - 2M 

This expression is zero if at least one of the following conditions hold: 

A = , Ur = , u = M/r . (30) 

The first condition {A = 0) can obviously be neglected because it implies that 
equation A3 disappears; the second condition means that u = a + bf{y), a, 
and b being arbitrary constants and f{y) an arbitrary function of the variable 
y. Nevertheless, that expression is not a solution of the system (121]), except 
for the case u = a, because Uy 7^ 0, and, since u = M/r has been forced to be 
a solution of A3, which is a non-linear equation, then the linear combination 
u = a + M/r is no longer a solution of A3. 

Hence, we must finally conclude that pr'^^W [A3] = iff u = M/r, i.e., 
whenever A^{x,u^"'^) = 0, as is the case, since that is the only solution of 
system 

□ 



Let {r, y} be an MSA system of coordinates such that (fT9l) is fulfilled for 
= 1. We can then state the following: 
Theorem 2 

The system of equations A^{x,u^'^^) = given by 

Ai(x,m(")) = Am = 

A2{X,U^''^)=Uyy = (31) 

A3(x,m(")) , 
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where Ai is the Laplace equation (with axial symmetry) (ITH]) . A2 the so-called 
supplementary equation [TD], and A3 equation admits a symmetry group 
whose infinitesimal generator is 

Odd 
or oy ou 



Proof: 

In [To] the null conditions on the prolongations of the vector field fl32|) act- 
ing on the first two equations of system (l3Ti) . pr^^V [Ai] = pr'^^W [A2] = 0, 
whenever these two equations are fulfilled, were satisfied. Now, we ex- 
plore whether pr*^^W [A3] vanishes for the solutions of the system (13T1) . It 
is straightforward to calculate (see [10] for details) that the second prolon- 
gation of vector (!32|) acting on equation A3 is 



pr 



(2)v [A. 



Urr (rAr + yAy - A) + Uyy {rBr + yBy - B) + 

Ury {2rCr + 2yCy - 2C) + 

Ur {rDr + yDy) + Uy {rEr + yEy) + 

^ {rCr + yCy) + 2^^^^^^2C 



2u 
2 



2u-l 
2 

'2u - 1 



(2u - I) 
(M. + ,A,) + 2^^--^A 
2m -2 

(rB„ + ,B,) + 2p— 



+ 



(33) 



Since pr*^^W [A3] = only needs to hold for the solutions of ( 13T|) . we can 
substitute the derivative w^r from equation Ai, 



2C 



D 



E 



Ur 



^^2u - 1) ^^"^ + ^'^y + 2CUrUy) - — - —Ur - jUy (34) 

into (133|) as follows 

pr^^W [A3] 



j/j, yvBr + yBy - B) + 
2rCr + 2yCy-2C (r^ + y^^ 



Ur 



A " A 

rDr + yDy + D{l-r^ + y^ 



+ 
+ 



14 



rEr + yEy + Ell-r^ + y^ 



2A 



(2m- 1)2_ 

(rBr + yBy) + 2- -B + r— ^ + 

2n - 1 ^ 2m - 1 2 ^ 2m - 1 V A ^ A y 



+ 
(35) 



By using expression in equations Ai and A2, the derivative u1 should 
satisfy the following equation 



2m- 1 ^ 



2 \Bul + 2Cm,m J C 



2 



Ur + 



A r2 



My . 

(36) 



By replacing expressions (l36l) . (!34l) and equation A2 in (l35ll we have that 
pr^^V [A3] = 

/ A. AA 2C 1 

+ 



M, 



2ra + 2?/C„ -2C\r^ + y^ 



A 



A 



Ur 



rDr + yDy + D{l~r^ + y^] + 
rEr + yEy + E{l-r^ + y^ 



2m - 1 
2A 



D 



r{2u - 1) 2m - 1 
2Ay E 



Uj'p'Uy 



4 



2m 



A J ' r2(2M- 1) 2m- 1 

AC f Ar A ' 

+ 



+ 
+ 



2B 



2u 



— (rB. + yB„) + — r-^ + ,-f 



A, 



(37) 



Since the coefficients of the various monomials in the first-order and 
second-order partial derivatives of u in this expression must be equal to zero, 
we impose the following conditions on the coefficients and their derivatives 
of equation A3 



C = 

A {rBr + yBy) = B (r A + yAy) 

2A 



D , , , , 2m - 2 
rDr + yDy-- (r A + yAy) + ^^^^ 



r(2M- 1) 
2Ay 



rEr ^yEy-\ (rAr + yAy) + = ^.^^^ _ . 



(38) 
(39) 

(40) 
(41) 
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The general solution of the system of equations Ai and A2 is as follows 

y C4 , . 

u = ci— + C2ry + + — . (42) 

Since {r, y} is a MSA system of coordinates (for the multipole order con- 
sidered) we can force equation A3 to possess the Monopole-Dipole Solution 
(among those from (H2|) ) represented by m = -M + \Mi (M and Mi be- 
ing the Monopole and Dipole moments respectively) and so, in addition to 
- (HTj) . we have the following condition {C = 0) 



A 



2 „ . 



2u 



-Bui + + Euy = {) . (43) 



Moreover, this imposed condition implies that the only solution of the system 
of equations (13T|) is the Monopole-Dipole m, and henceforth we can replace u 
by u in equations (155]) - (HT]) . leading to the following expressions: 

(a) + yj-'y = 

2m - 2 9 

[0) TKj. + yKy H K 



2u - 1 r{2u - 1) 

2m -2 2y 

where u = ^, k = ^ and /i = ^. Equation (jlUa) requires u to be an 
arbitrary function oi y/r and the solutions of equations (jBlb), (jlUc) for k, 
/i are the following functions 

.-u(v/r) ^_JZ±IMA 2y-F^ 

Equivalently, these solutions of the determining equations ( 14^ can be written 
as follows: 

vA = B 
(2r + Fi)A = {r^ -2Mr -2Miy)D 
{2y - F2)A = -((r^ - 2Mr - 2Miy)E . (46) 

These functions z/, k and /x are related by means of condition (1431) . leading 
to the following relation between Fi, F2 and v 

M1F2 = 2M2 + 4:M^ (-] ^ + 8MM1 - - 2MiV + MFi + 2MiFi - . (47) 

Vr / r r 
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The first equation of (H^ and equation (H7I) allow us to write the following 
expression: 



{2y-F2)A = 2MiB+A 



y ( 

Ml r \ 



r) \Mi r 



(48) 

Finally, by substituting expression (HHll . and the function Fi obtained from 
the second equation of (H6l) in the third equation of (H6l) . we have that: 



A 



..y ... fVV M 

6v-2 8M--4Mi - +2r — 

^ Ml r \rJ Ml 



+ 52Mi + 



-D 



-2 2Mr - 2Mi2/) f ^ + 2^ 
\Mi r 



+ E [(r^ - 2Mr - 2Mi?/)] = . 

(49) 



We must therefore say that the second prolongation of the vector v acting 
on A3 vanishes whenever the system of equations (13T]) is fulfilled, iff condition 
(149|) holds, and hence the proof of this theorem can be concluded since this 
condition poj) is equivalenio to saying that the Monopole-Dipole function u is 
a solution of A3; that is, the assumption from the beginning of the theorem 
if coordinates used {r, y} are MSA coordinates. 

□ 



The relevance of the theorem comes from the relationship that can be 
established between the existence of the symmetry and the gauge of coordi- 
nates that provides the Newtonian form of the Monopole-Dipole Solution. 

If we recall definitions fllOtllip . then the last expression (H9!) . in addition 
to the first equation (14T|) (C = 0), provides the following explicit conditions 
in the coordinate transformation: 



(a) LBi{r, y) = R^rnyn + (1 - uj^)r^y^ = 

(6) 2LBi{r, r) [-Mr^{r - M) + Miy{AMr - Sr^ + 2yMi)\ = 

r{r^ - 2Mr - 2Miy) MirLB2{y) - (Mr + 2yMi)LB2{r) 

2MyLBi{y,y) . 



+ 



(50) 



Note that equation ([43|) is exactly equal to expression (j49|l . 
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Note that if we take Mi = in the above expression, we obtain condition 
f l27p . which must be used for the determination of the coordinates in the 
Monopole case. We shall discuss these results in the following section. 

3.2 Characterization of the MSA systems of coordi- 
nates 

A) The Monopole Solution 

From theorem 1, one can conclude that there exists a symmetry of the 
system of differential equations (12T1) iff we force equation (!23|) to posses a 
solution of the Monopole type u = M/r; in other words, equation (!27|) must 
be satisfied. Equation (l23l) is the Ernst equation for the function u with the 
constraints given by the other equations of the system (1^ . and it should 
be taken into account that the system of MSA coordinates that we are using 
allows us to characterize the relativistic Monopole Solution with a function 
u written as the Newtonian Monopole. Therefore, since the only solution 
of system fl2T]) is the Monopole Solution, we can state that this system of 
differential equations admits a symmetry group that can be related to the 
uniqueness of the solution of the system by means of the existence of MSA 
coordinates. 

Furthermore, condition (127|) allows us to determine the gauge of coordi- 
nates in which the relativistic solution having only the Monopole Moment is 
given by a function with the same analytic form as the classical Monopole 
potential in NG. We proceed to do this in the following way. First, we substi- 
tute the coordinate transformation (1121) in condition fl27p . which becomes a 
constraint over the coordinate transformation, taking into account definitions 
(HOKID; i.e. 

r(r - 2M)LB2{r) = 2(r - M)LBi{r, r) . (51) 

Second, we solve the corresponding differential equations that appear at each 
order in the power series expansion, and the uniqueness of the solution is 
provided by the following boundary conditions: 

r{uj = ±l)=R + M, (52) 

which implies that all functions fn{^^) vanish for all n > 1 along the axis 
orthogonal to the equatorial plane. 
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As already noted, that radial coordinate r is merely the Schwarzschild 
standard coordinate, the coordinate y being free of constraints because of 
the spherical symmetry. The system of coordinates characterized by solving 
equation flSTl) with boundary conditions fl52|) . can aptly be said to be adapted 
to the Monopole symmetry group, whose infinitesimal generator is fl20|) . for 
several reasons. First, we see that, written in these coordinates, the solution 
does not depend on the angular coordinate. However, in addition another 
feature contributes to characterizing these coordinates: we refer to it as 
MSA because of the interrelation between the existence of the symmetry 
and the system of coordinates itself. Second, the function u that describes 
the relativistic solution with a finite number of RMM acquires the form of 
the classical Multipole potential, and hence all the conclusions obtained for 
the Newtonian case can be assumed again for this function, which can be 
considered as the group- invariant solution of a system of differential equations 
f[T8|) that admits the symmetry. 

B)The Monopole-Dipole symmetry 

In analogy with the previous case. Theorem 2 allows us to establish a 
relationship between the existence of a symmetry of a certain system of 
equations and the MSA coordinates for the Monopole-Dipole Solution. We 
have seen that the uniqueness of the solution of the system of differential 
equations fl3T|) can be deduced if the Ernst equation for the function u written 
in MSA coordinates is required to have a solution with the analytic form of 
the classical Monopole-Dipole gravitational potential. At the same time, 
this condition leads to the existence of a symmetry group for that system of 
equations. 

The main goal obtained in the previous case (Monopole) is the calculation 
of the coordinate r by means of the constraint flSlI) that Theorem 1 intro- 
duces in the transformation of the coordinates; with appropriate boundary 
conditions fl52|) . the choice of the functions /n(cu) is unique and r is fully 
determined. 

Nevertheless, we cannot perform the complete determination of the MSA 
coordinates for the Monopole-Dipole case by using the constraints fISU]) in- 
troduced by Theorem 2. Equation fISPl a) means that the new coordinates r 
and y preserve the orthogonality since they must be asymptotically Cartesian 
coordinates. This equation also implies that there are not cross terms in the 
metric written in MSA coordinates {gij = 0, for i ^ f). If we substitute the 
prescribed gauge transformation from Weyl coordinates fll2l) in constraint 
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(pTla). the corresponding series expansion leads to the following equations 

k fc-1 

(1 - uj^) f-{uj)g'k^^{uj) = kgk{uj) ^(n - l){k - n)/„(^)^fc_„(^) , (53) 

i=l n=2 

where goioj) = 0, the symbol (') denotes the derivative with respect to the 
variable u, and k is the different order of the expansion in the parameter 1/R. 
This expression allows us to write any function gk{oj) in terms of the functions 
fn{i^) and gn{'^) of lower order {n < k) and their derivatives. As can be seen, 
the good behaviour of the functions gni,^) is recovered, since they must be 
zero along the axis of symmetry (a; = ±1). With the expressions of gk{oj) 
obtained from (1531) . we may solve the corresponding equations at each order 
of the series expansion of the other constraint, (ISUlb). for the functions fniyj) 
alone. But now we do not have a suitable boundary condition to obtain a 
unique solution of the functions fnil^) for the complete determination of the 
coordinate r for this Monopole-Dipole case. We can demand that the limit 
Ml = must lead to the same functions as the Monopole case, but this only 
allows us to determine the arbitrary constants of integration for /i(w) and 
/2(c(j), because from the next order onwards the equations involve arbitrary 
constants, since the functions /n(co') for odd n are null in the Monopole case. 

With this procedure we obtain a family of coordinates that transforms 
the Ernst equation into another one that admits the function m, representing 
the Monopole-Dipole Solution, as a solution. Nevertheless, the coordinate 
transformation is not unique. 

4 Conclusion 

The exterior gravitational field of an isolated and static compact body with 
axial symmetry is described in GR by means of the Weyl family of solutions, 
which depends on a set of arbitrary coefficients {a„} whose values univocally 
determine each specific solution. If one is looking for solutions that are well 
known and physically meaningful, it is necessary to relate this set of coeffi- 
cients to the RMM in order to make a suitable selection of them. If we work 
with an MSA system of coordinates, then the solution that describes the 
gravitational behaviour of compact bodies with a prescribed multipole struc- 
ture can be constructed by identifying the function u with the Newtonian 
potential and considering the constants of the classical potential to be ex- 
actly the RMM of the solution. The transformation of a solution (Weyl) into 
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another one (MSA) fixes the change of coordinates by requiring the function 
M to be a solution of the corresponding Ernst equation. 

Nevertheless, this procedure needs, a priori, to know the set of coefficients 
{an} of the desired Multipole Solution, although the existence and explicit 
knowledge of these systems of coordinates is relevant enough and they be- 
come very useful, at least for the following topics: Application of this work 
could shed light on study of the influence and relevance of different RMM 
in the behaviour of test particles along geodesies [16] for different sources, 
since the MSA coordinates provide us with the exact Multipole Solutions 
in a very simple way. The deviation of the source from the spherical con- 
figuration is a very important feature, for example, for describing the fate 
of the collapse of self-gravitating systems [T7]. The calculation of circular 
geodesies at successive distances from the source can be used to determine 
its multipole structure. Additionally, some authors have attempted to [I5] 
relate the RMM to the structure of the source by means of quantities defined 
over the distribution tensor of the source. The existence of MSA coordinates 
for any Multipole Solution seems to be a very useful tool for achieving these 
aims in the frame of global stationary axisymmetric solutions of the Einstein 
equations. 

Except for the Monopole case, in which the gauge is well-known (Schwarzs- 
child radial coordinate), the MSA system of coordinates for the other cases 
are given by means of two series expansions in the inverse radial Weyl coor- 
dinate. However, if we work at large distances from the source, the approxi- 
mate character of the coordinates is negligible and the change of coordinates 
is completely determined. In addition, the function u is an exact solution 
with the finite number of desired RMM. Also, we have defined a family of 
static and axially symmetric exact vacuum solutions with a prescribed mul- 
tipole structure in a system of coordinates defined with a suitable order of 
approximation. 

Finally, this work affords another conclusion: two theorems have been 
proved that allow us to establish a relationship between the existence of 
a certain symmetry of a system of differential equations (the correponding 
Ernst equation for the function u is included among them) and the existence 
of a system of coordinates in which that function u can be written analytically 
equal to the Newtonian potential but in terms of the RMM. 

This result is relevant in itself, and some implications can be derived from 
it. In particular, the construction of the MSA system of coordinates for the 
Monopole case is supported by the proof of Theorem 1, without knowledge of 
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the corresponding set of coefficients {a^}, since the existence of the function 
u, or the MSA system of coordinates, is equivalent to the satisfaction of the 
corresponding Ernst equation by that function. This condition requires that 
the functions involved in the change of coordinates must satisfy some differ- 
ential equations whose solution is unique for suitable boundary conditions. 
Accordingly, the existence of a one-parameter group of transformations can 
be stated, whose infinitesimal generator is fl2U]) . which represents a symmetry 
of the system of equations joined by the Laplace equation with axial sym- 
metry, the supplementary equation and the Ernst equation for the function 
u written in a system of coordinates adapted to that symmetry (MSA). 

For the Monopole-Dipole case. Theorem 2 allow us to establish the same 
relationship between the existence of the Monopole-Dipole symmetry and the 
corresponding MSA coordinate system, although unfortunately the charac- 
terization of the gauge by means of the conditions provided by the theorem 
leads to a family of undefined coordinates in terms of arbitrary constants. 



5 Appendix A 



The following expressions show the first ten RMM of any Weyl solution in 
terms of its coefficients a„: 

Mo = -ao , Ml = -ai , M2 = -02 + , M3 = -as + aial 
8 2 19 5 6 2 

Ma = -aoOn On H CtnOi — Ua 

7 ° 105 ° 7 ^ 
4 19 12 2 

M5 = -a-sal - — aiOo + ya2aoai - cts + ^^i 

20 34 o 2 23 4 389 7 17 2 60 2 

6 2 
+ ^020^ - ae 

206 o 2 389 120 2 595 . 126 2 23 2 

^ 143 ^ ° 495 ° 143 ^ 429 ° 143 ^ 13 ° 

282 20 1504 3 

+ Y^^^ctoC^i — a? + —o.3aQa2 — ^^020.10.0 

2948 2 ^ 652 , 40 o 257 „ 226 

Ms = aiOn Ctn^i H 0-9 On H 040002 + 

1365 ^ ° 1001 ^ 143 ^ 3465 ° 143 

4464 2 2 1744 3 5204 2 3 28 2 58 4 
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+ 



+ 



12 2 2 44312 . 240 28 

—040^ + 20600 + 43Q4g Q2ao + ^030201 + jgOsaoOi - Os 

163508 4 , 257 o 3192 . 1988 , 120 , 
^^aoo, - — a,ao + ^^a,a, + ^^a,a^ + —a,a, + 

486 4 216 2 11724 3 40 44424 3 

^^^5«o + ^^^-'sO'i - 2431 '^'^'^^'^0 ~^ Y^«6«o«i - ]^7Q]^7'^2ao«i + 

10908 5 3430 11944 . 2 9080 3 366 

85085^^ + 2431"^"°"^ " W"^"^"-^ " 2431^^"^"° + 22l"^"'^"^ + 
418144 5 38 2 378 5748 2 2 

85085 '^^'^I'^o + Yy«7«o + ^^040201 - 09 - Y309'^2'^i'^o 

828 10041124 , 4 10908 4 17389 o 

10 ^" ^ 323 969969 ^ ° 17017 ^ 20349 ° 

555820 2 7 226580 2 5 459700 3 2 262556 2 3 

223839''i''° + W79"''2"° ~ 323323''2''° " 13^ 

3500 2 12902 f. 39150 2 30870 2 2624 , 

+ 4199"^"^ + 7293""^"° + 46189^^"^ + 46189"°"^ " + 
336 2 47 2 193130 3 75180 896 3 

+ 323"^"^ + 19"«"° - 46189""^"^"° + 46189"^"^"^ " 323^^"'^"^ + 
94632 2 2 70412 , . 566 257460 2 2 

- 24871"^"°"^ + 24871"^"° + 323^^"°"^ " 16189-"^"^"° + 
5992 427568 2 44152 , 7416 

050003 03020iOn H OsOiOn H O5O2O1 + 

4199 46189 ^ ' ^ 7293 ° 4199 

443699 y, 24116 o 

+ 8729721^° - 119^"^"^"° ^'^^ 

From the above expressions we can extract, at each order, the correspond- 
ing coefficient o„ in terms of the RMM: 

oo = -Mo , oi = -Ml , 02 = —Ml - M2 , 03 = -MiM^ - M3 

8 16 
04 = --Mo^M2 - -M^ - -MoM^^ - M4 

as = -^M2M3 - ^MoMiM2 - MiMo^ - ^Mf - Ms 

Of I 

06 = -^MiMoMg - ^Mo^M2 - ^M^M^ - iZMo2M4 - ^MoM| + 

- -MlM^--Ml-M^ 

fi 128 n 49 4 120 , 18 . n 

07 = -MiMo^ - — M3MiM2 - -Mo^M3 - — M|Mi - -MfM^ + 



+ 
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23 o 126 o 282 20 

—M^M. MfMs M0M1M4 - My M0M3M2 

13 ° 143 ^ 143 7 3 2 

, 28 , 820 , o 92 , . 108 . 
-2Mo2M6 - 39MIM0 - 429^2^0 - 33^1 ^0 - 143^°^' + 

Im^ - — - —M^M^M^ - l^MsM^Mi - — Mn^M4 + 
9 ° 33 ° 143 ° ^ 429 ° 143 ° 

40 o 226 28 240 

Mg - ^43^2 - Y^MoMiM2 - — M0M1M5 - — M2M1M3 + 

12 , 

13 ^ 

MIM3 - Mg M0M5M2 Mf MfM^Ms + 

143 " ^ 221 ° ^ ^ 715 ^ 143 ^ ° ^ 

432 o 888 , 224 . 30 , 

143 ^ 143 ° 143 ° 13 ° 

— M5 - —M^MiM^ - — M2M1M4 - —M^M^Mi + 
221 ^ ^ 143 ° 221 143 ° ^ 

38 , ,2 . . 40 , , , , , , 604 , ,4 , o 3430 ,,,,,, 
—M^Mr - -MoMiMe - r^M^M^ - ^M^M^M^ + 

143 ° ^ 2431 ^ ° 
336 , ,2 , , 14940 , ,2 , ,2 , . 11010 ,,,,,,, 
323 ^ 2431 ^ ° 2431 « ^ 

^^Mn^MiMs - '^^"^^^ MiM4M3 - ^^MoMsMg - — MnMe + 
221 ^ ^ 46189 4199 " ^ 17 ° 

7416 294 f. 5012 o , 47 , 

^M,M5M2 - — Mo^M4 - ^M3M3^ - -MlM, + 

tS^'^° - |^AfoM2M6 - I^MoM.M, - ^M,^M2 + 

10728 , ,2,, ,,2 ,, 540, ,2, ,7 30870 ,, ,,2 

M^MqMo - Mio M?MI MqMI + 

2431 1 " 2 10 1 4g^gg 4 -r 

460 , , ,2 176 , , , ,2 , . , . 39150 , ,2 , . 

143 ° ^ 17 ° 46189 ^ 

148,^2,^4,^ 1064,^,.,^ 7728 ,^3,^,^ 

^^102 143 1 3 2431 1 3 ^ 

175,,, „ 300 ,,2,,s 1 ,.11 3500 , ,,,2 

M2M^ M^MS M}^ M2MI (55) 

143 ° 187 ° ^ 11 ° 4199 ^ ^ ^ 
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6 Appendix B 

The following expressions show the functions fn{'^) appearing in (fT2|) . for 
the Monopole case. Since the spherical symmetry only requires to define the 
radial coordinate, and the function u does not depend on y, the condition 
(EO) provides itself this result: 

/iH = M, f^[u) = -]^M\~l+u^) , /3H = 

Uuo) = -— M6(-35cj^ - 1 + Ibu^ + 21u^) , fyiu) = 
16 

/8(cj) = _^M^(5 + 429u;S - 924^;^ + 630a;^ - UOu^) , fg{uj) = 
128 

/lo(cu) = — M^°(7 - 2431cu^° + 6435cj^ - 6006cj6 + 2310cj^ - 315^^) 
256 

(56) 

For a more general case, we have calculated the MSA coordinates for the 
solution having only the Monopole, Dipole and Cuadrupole moments, and 
the functions /„(cc;) and 5'„(a;), up to order 10, are the following: 



f\{u) = Mo , /2M = -^Mo^ ^ 1^2 

1 

/3M = MoMiu{l - u^) + -M2{1 - 3uj^) + -^tu^ 



2 ' 'Mr 







5 , .2 2 8 M^ 5 _ , . 39 2 , . , . M: 



/4H = -^MV- 3^^'-4MoM2^^ + ^a;^MoM2-3Mia;^ + 

+ 1m,^ + iMo^(6..^ - 1 - + TMic^^ - ^M^^ + 
19 , , , , 2 Mi^ 
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1 103 _2 4 M2 87 _ _ o 31 _2 2 ^2 

85 9 , 

+ — M1M2UJ MiM20J^ 

28 4 

N 13 3293 ,^2 2,. 2 1777,, 2,. 2 4 223 , ,2 , .2 

■'^^ ^ 105 840 ^ ° 168 ° ^ 840 ° ^ 

21 ,.6 6 46 ,^ Mf 1247, . 1 , . 

16 ° 5 Ml 616 2 16 ° 

308 ,Mf 1378 . 63 _2 « 214 Mfa;^ 



26Mfu; 15 _2 fi 13089 _2 4 10347 2 

^ M^uj^ M^uj^ H + 

15 Mn^ 8 ^ 616 ^ 616 ^ 



'0 



45,^3,, 6 .,.3 5 496 MV 25 MfMs 

- —M^M2UJ^ - AM^uj^ H ^- H i — ^ + 

8 ° ^ 105 28 Mo 

6 , . , . , . 901 M2MV 713 M2MV 5 

- -Mi^MoM2 + — ^ttI — +15MoMiM2^^ + 

5 10 10 

6597M2u;4M2 247 M| Mi 366 M| Mi cj^ 
^ TIT Mo 20 M2 ^ m| 

1617 M|Mia;5 1781 MVM2 35 2 3 

" ^0 Ml 70 M^+16'^^°"15^^"' + 

+ Am^o; - f^u^M', - ^-^MlM^u' + ^M3M2a;^ + ^M3M2 
29 

- — Miw^MoMs 
5 

142687 17 
f,{uj) = Mi^c^^Mo^ - -^^Mfc^^Mo + —M^M^u^ - 2M^M2u;^ + 

I006O z 



98 MV 475 Mfa;^ UA Mlcu^ 8503 Mf 



^ 9 M5 9 M5 5 M2 315 M^ ^ 

33,,,, 3 7 99,.,^ 7 2618Mfo^6 88049 ,^ 3,, 2,. 
- -MoM.V - -Mo^Mio;^ + _ ___mVMo^M2 + 

+ tSS^^^i^^o - ^^0^1^' - ^Mo^Af2a;' + bM^M^uj'^ + 
iooOU o Z 
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839 Mfu^ 49Mfuj^ 107 MJ^Ma 191 MfM^ 

630 Mo Mo ^60 60 M^ ^ 

1889 , 107Mfcu 8783M|cu^M2 15329 M^ 

414 M^iu^ 123 p 73 Mi-^cj^ 981 M^u^ 

^ M?M2 H — + ^ + 

5 Mo^ 140 ^ 126 Mo 10 Mo^ 

189,^5,^ 5 81^,^2 6 3 M| 2 Mf 

+ M^Micu^ H MqMIu^ H 1 H + 

8 ° 2 " ^ 2 M2 15 Mo 

13 Mf* 61 , 9 4 q 1583 , , 

56 ° 60 M^ 4 M3 

3947 MaMi^u;^ 897742 Miu;3^| 4147313 Miu;5M| 
~60 M| 5005 Mo ^ 20020 Mq ^ 

16793 Ml cjM| 209 ,^2,.,. 7 121,^2,. 6 
+ Mo - ^^o^iM^c.^ - — M,-M2C.^ + 

1303 M|MV TGSlMaMfo;^ 2155 MsM^ 
^ ~V1 M| 35 M| mI ^ 

15 , , , 2711 ,,,,2 4 38 , , , ,2 2 159 MaMfo; 

+ —M^MiUj M0M2V MqMW 

8 ° 140 " 5 ^ 10 M2 

23719 _ _2,. 2047 _2 a.. 18113 MaMfu;^ 
9240 ° 420 ^ 60 

/3(^) ^ ^^2^8 l^M«u;« - — M^a;^ + — M«u;6 + 
y 32 -^28 ° 64 ° 32 ° 

1318143 , ,4 4 26323 , « 927191 , 2 383 , , , 

112112 ^ 840 ^ 840840 ^ 1056 ° 

138851 M| 127 Mf 48658 Mfw^ 227 MfM2 

24024 Mo ^ 315 M^ ~ 2205 ~ 70 M^ ^ 

936731 , ,2 , .2 3655 , , 19486 1156 Mjto^ 

- Mn^M,^ M^Mf + — ^ + 

672672 ° ^ 7392 ° ^ 105 Mq 21 Mq^ 

I66M70; _ 117^4 8_ 13553439 M|o;^ 62329 MfMf 

^ 45 M6 ~ "32" 40040 MfT^ 8008 M^ ^ 

431143 333 

+ M0M2M? M^u;^MiM2 + 27M^u;'^MiM2 + 

1681680 ° 2 1 35 ^ 2^ ^ 2^ 



27 



- —MlJ'M^M^ + —LuM^MiM^ + -^Mluj^'MoM^ 

52341 M|a;2 892359 M|a;6 292162 Mfa;^ Mfa;^ 

^ ^20 Mo~^ 3080 2205 ~ ^ 

1496 MV 27494 MV 3662 Mfu; 4897,,, 

^ ^ M^uj^M2 + 
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+ ^MoVm2 + + y^m^^'m! - —M^cu'm! + 

_ ^^MoVM^ - ™Mo^cu8^2 6^2 

560 ° ^ 32 ° 1 ^ 84 ° ^ 
10 ,.2,.s 15916 2,. 4,. 2 85223 , 5 , 

21 ° ^ 1155 ° ^ 1121120 ^ 128 ° 



385723 Mf g;^ 34096M^ _ 715 ^ 
2205 + 2205 " ^^o^W^o; + 
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+ H M^u^M^ H ^-^^ h 

735 Mo 980 ^ 35 Mq^ 

29103 Mfuj^M2 56453 MfMiCj^ 38004 M^M^uj^ 

~245 Mo ^^40 M| 35 Mf + 

6221 ,, ,,,2 1077,, 7,, 2 663 ,,,,,, 2 8 

+ Miu^M^ Miu^M^ MoMsMf u;^ + 

245 ^ 4 2 16 ^ 

22771 M|MV 34285679 4 , , ,2 1143841 . , .2 w 

30 M^ 560560 ° ^ 140140 ° " 

112353 M|cu^Mi 35619 Mscu^Mf 4121 6,^2,^2 

-I- 7: h = — OJ Mn Mo ~l~ 

140 M3 35 Ml 77 02^ 

5858 ujM^Mf 904 UJM2MI 1123 uj^M2Mf 

+ 105 M^ 35 Ml + ^ m| + 

374531797 M|M2(^4 M^uj'^ M2 73619837 MlM^a;^ 

_(_ _(_ 265 - + 

280280 M2 Mq 56056 M^ 



2 



28 



695879 M2Mfu^ 125A99 MWMf 18163 _ , . 

+ t^^mT - ^''^^^'^ 

^, , 3786191 ^3 3,^3 105,^ 3,^7 693 ,^ .,^7 
/qH = 72072 + — ^i^ ^o -^^i^X + 

+ 9M2cc;X' + — M,VMo'M2 + 3655803 4 3 
1 ° ^ 40040 ° 

2139, . 2, .2,. 848053,^,^4,^ 1714497,^,^4,^ 5 

- M^uJ^Ml M2 Ml Mn Mao; Mi Mf Mao;^ + 

140 1 02 ggQgg 1 2 g^gQ 1 2 T 

401,, 2,. s 2 1800157 , ,0 , ,3 , 354703,^3,, 3 

20 ^ ° 9240 ° ^ 90090 ° ^ 

35 1 

+ 2M2M0V - 16M2M^uj^ MiwMj + -MIMIM2 + 

16 7 

253 4 2,. 15772 Mfcu8^2 752747631 M|cj^Mi 
+ — MicMo-^^ ^ + 

14495 _3_3 7 10449 M|a;2 22469 M|M? 



1306727 Mi^a;5 13149 M|a;^ 4301 Mf, 



2205 28 M^ 9 Mj ^ ° 

934429 Macu^Mf 1271,^4 4,, 136,. 2 ,^3 

—ir-^ MVMo + Ml MlMn^ + 

280 M6 315 ^ ° 11 2 21 ^ ° 

231,, 3,, 2 8 30763 Mfa; « 119,. ,.4 « 

+ M^Mlu^ + ^— + 34Mn^M2a;^ H MoM>^ + 

2 ° ^ 84084 Mo ° 3 ^ 

1807363 M/cj^ 113867 Mfa;^ 20725 Mfa;^ 

^ ^ 2940 M| 180 M^^~84 M^ ^ 

608MfM2 1621 M|Mi^ 5189 MfM^ 35887 M^ 
^ 147 Ml ^ 168 M^ 560 M^ 1260 Mj ^ 

210779 Mfo;^ 50191 M^ 71303 M/o; 769 MfMa 

^ 1260 M^ ^ 2520 Mq^ ~ 8820 M^ ~ 210 M^ ^ 

20085 , ,4 , , , , 7 5805579 MWmI 662141 Mjcu=^ 

+ Mn^MiMaCJ^ ^ , ^ + ^ + 

56 ° 1120 M^ 4410 Mq^ 



819555 M^u^Mf 578719153 M|Mf 



a; 



+ ^ , ^ 105Mn^M2M>** + 

112 M^ 3363360 M(f 2 1 ^ 

6673937 M2Mfa;3 8809390459 M^M^ 405 , ,3,, „ 
5880 Mq^ 3363360 M^ 16 ^ ^ 
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9431951569 MlMfcu^ 12523453 _ 429 ,,7 _ 7 
1121120 M3 1681680 ^ ^ ^ 4 ° 



715,, 7,, 9 316,^,. 6 1011517 Mfu; 

16 ° 315 ^ 8820 Ml 

7761,^3 8 845 , 9 85089 , . g 4086767, . . 

20 ^ 8 ° ^ 2695 ^ 10780 ^ 

42251 MfcjS 1148139 Mfcu^ 1125,,,,,, 2 q 

+ Kr- + MqMiM^uo^ + 

504 140140 Mo 16 ^ 

18563603 Mluj^ 2816 Mfa;^ 132543 M^u^ 

252252 MfT ^ 140 M|~ 

20493 M|a;8 ^.,,,6,. 8 437691 ^ „ ..,,5,, 2 8 
+ ^5 " 21Mo^M2a;« - ^^Mfo;^ + MMq^MV + 

, , 2546307 o , 2536 , , 

79123 3 7 11393363 :172^'Mf 1259203 ilZ/^M/s 

" ^0^^^'^^^"^940 M^ + ^470 + 

990122241 AffMicu^ 33193 . . ^ 259 , , 

+ MIM^uo^M2 M^Mlu/ + 

224224 M2 420 ° ^ 2^20^ 

5655 , g 35031 Macu^Mf 69329927 M|Micj 

- M^MiM2U^ ^ ^ H ^ / + 

32 ° 98 M2 480480 

126526 MaMi^cu^ 12186841 MsMfcu^ 1353 „ . 

— + M^MW + 

735 2940 28 ° ^ 

164152063 M|Mic^3 2991113 MlMfa;^ 3097M|M2a;« 

101920 M| ^ 10780 Mo Mo ^ 

42701053 , ^ , ^ , . 19308785 MfM^o;^ 

+ MoMicuMi ^— h 

1681680 ^ 12936 Mq 

951347 MlMi^u;^ 388861 MsMftu 186671 MfMaa;^ 
1680 Ml ^ 5880 M^ ^ 840 M| ^ 
31 4 341 Mf 9Af| 314309 MfMacj^ 

MnMn TT -\- TT T\ -|- 

315 " ^ 735 M3 4M3 168 M^ 

83 Mf 4065127 Mlcu^M^ 40869853,,,, ,,,0 

+ \ — 7-^ H MoMia;^Mp^ + 

180 Mj 1120 40040 ^ 

8890071 M^uj^Ml 5957607 M^uj^Mf 653883 , , , ^ , , . 

+ 1120 Ml - 560 Mo^ +^^"^^"^- + 



4 
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557573 MlM^u^ 7630153 MlMfuj^ 28687558 _ _ . . 

+ — + MqMiuj-'MI + 

1120 1680 35035 ^ 

^ 3497168 MlMfcu^ ^ 722625097 M^uj^Mf ^ 1294039 MsMfa;^ 



8085 Mo 101920 280 

, , , 65998679 _ , . 315 , 2 1522198 Mf tu^ 
/loH = ^ 001 1 ono ^oM| - — Mo^V - ^ 



3811808 ^ 256 ^ 2205 Mq^ 

1091827871 MlMiuj^ 116182867 MfM2Uj^ 

168168 Mo ^ 1261260 Mq ^ 
823063229 . „ 1071 , ,2 , .2 in ^.^,,.2,.2,. q 
140140 " 16 ^ ^ 021 

28295301319 . . . 86424015889 M^cu^Mf 

+ Ml Mno;^ H ^ — ^ + 

6126120 ^ ^ 3363360 M^ 

90647147471 ^2,^2 2 41574419,^2,^2 4 

257297040 ^ ^ 2144142 ^ ^ 
510812297 Mfu;6^2 720726269 Mfa;M| 
26460 Ml ^ 2522520 M^ ^ 

24253205017 Mftu^Ml 67749616547 M^cj^^l 
^ 1441440 Ml 30270240 M| ^ 

37401809 MVM2 57147915281 Mfu^M^ 
^ 5292 M| 3027024 M^ ^ 

8637393 M^uj^M^ 4462345 1 M^uj'^ M2 

1120 M^ ^ 17640 Mf ^ 
2124423 M^uj^Ml 165 9717 M^u^M^ 
^ 160 M^ 224 M| ^ 

16666 153 -l/j^^'^-l/ M2 1493281 1329 

10584 M| ^ 315315 M^ 
1482330167 Mfuo'^M'^ 201899419 Mf w^Mg 
^ 687960 mJ ^ 13230 M| ^ 
47883811 Mf cu^Ma 49533 MicuM| 

8820 Ml 560 M^ ^ 

25556833901 Mfuj^Ml 240257 Mf cu^Ma 

672672 M| "^~60 Mf ^ 
1801203 MiUj^M^ 2995594751 Mfu^M^ 
^ 1120 mJ ^ 252252 M| ^ 



+ 
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305701 Mfuo^M^ 407164048 Mfu'^M^ 
^0 Ml 105105 Ml ^ 

18105, in 23308919 Mftu^Ml 

- M^M^M2UJ^ H — ^ + 

32 ° ^ 840 

11923937 Mfcu^Ml 13639837 M/cu^Mg 

420 Ml ^ 1260 Ml ^ 

2442373 M5a;M2 _ 983324 MV^2 

26460 Ml 1323 M| ^ 

5623987 Mlu^M2 6494269 Mjc^^Ms 15103 Mf M2 

3780 iW^^ ^ 756 iW^^ ^ 1470 M^ ^ 

33154333 Mfa;3M| 209051 Mi^a;M| 29207791 ^2 2,.s,. 

+ —ir^ —ir^ Mfuj^M^M2 + 

7560 M^ 1080 M^ 550368 ^ ° 

235639 M-fcuM| 12025 Mi^tuMs 8260429 , _n _n _ 

+ + uj^M^M^Mi + 

1008 Ml 189 Mq^ 32340 021^ 

65076171 Mo^a;^ 56525 , 2664951647,^0 6..-^.. 

+ H Mn^Mi Mtuj^M^M2 + 

160160 M2 82368 ° ^ 1681680 ^ ° 

465748042,^2 4,. 3,, 3401 ,^ 3,. 5,. 259 ,^ 7,, 5,. 
+ n.^n.^ MV^o^2 + -—M^u^MlM^ + —M^u'M^M^ + 
y4o94o iOo o 

3296 ,^ 156327 M|Mia;9 982936967 MfM^ 

- MiUj^MlM2 + ^ — i + 

35 ° 40 Mo 30270240 Mq^ 

29351 „ 8, 42882109,0 5,^,^ 

+ ——Mlu^MlM2 - ———Mlu^MoM2 + 
16 iZJbb 

4225037,^3 3,.,. 604591 ,^3 7,^,^ 

+ Mfa;^MoM2 H Mfa;^MoM2 + 

48510 ^ 1470 ^ 

50027 ,,2,. 2,. 24439994 , , ,2 , .2 , . 21011 ,,,,,,3 

- cjMn^M2^Mi + u^M^M^Mi o^MoMaMf + 

9702 ° 2 24255 ° ^ 72765 ^ 

212 479261 

+ 28MiUJ^MlM2 M^MiM2UJ cu^M^M^Mi + 

° 105 ° 980 021 

7 81941597 Mi^Mau;^ 128708191 M|Mia;^ 

+ MA^ — \ — h 

256 ° 1513512 Mq 105105 Mq 

60156797 1452013 ,^5 3 290663 ,^5 , 

181621440 ° ^ 72765 ^ 2205 ^ 

2346121 , . 5 44089, . 546 , . . 

+ ^425^^^^+ 72765^^^-^^^^^°^^ 

+ ^MVM,^-^M3.3M,^-^MVMo^ + 
105 ^ ° 315 ^ ° 1120 ^ ° 
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106269 ,^4 2246545753 ,^4 2.. 2 799877851 ,4 fi,.2 

448 ^ ° 181621440 ^ ° 3363360 ^ ° 
1744579433, . 4, . 8444941 „ 25144057 4,, 4,, 2 

18162144 ^ ° 36960 ^ ° 1681680 ° ^ 

893555 ,^2 R,.fi 1422352 n 21271059 ,0 « 

1344 ^ ° 35035 ° ^ 12320 ^ 

12075043 M|M2 1275,., „ in 
308880 Ml 8 ^ ^ 

65870521 Mfa;6 352 , . « 11235935 

+ Mfuj^ ^- + 

26460 3 ^ 5292 

876456171^1^0;^ 5873111 M^^Ma 27874 M/a;^ 



19845 1681680 Mq 45 Mq^ 

38337021 M^uj^ 10355481 M^u;^ 13368209 



+ 



12320 Af2 ^ 2080 840840 ° ^ 

17022571101 ^ g 780547 4,.,^ 

9529520 ° ^ 7392 ° ^ 

214, ,4, ,3 112 o,,4 57863 g,.?,. 
315 ° ^ ?>Ml ° 224 ° 
2866051 2,. 7,. 174555 8,^4,, 2 13022941 ^ 

192192 ° 2464 ° ^ 15135120 01^^ 

88845511 , . 2,. 6 114075 8,^7,, 18619, 

2882880 ^ ° 448 02^ 20384 ° ^ 
2935 ,^,^7 8415,, 2,. 4 in 5525 ,,.,, ^ 
9152 ° 64 ° ^ 64 ° 

16575, .,,2 in 767417513 , ,2 , .2 78404639 

- MXM?u^^ H Ml Mo ^- + 

64 ° ^ 64324260 ^ ^ 3783780 M^ 

6585608633 , , , ,0 2 929535058 Mf o^^ 

57177120 " ^ 2837835 M^ 

1275 ,,,,, in 18322867 M/o;^ 189551 M^ 

- MoM|u;^° H ^ + 

32 " ^ 26460 34398 M^ 

2690985929 ,,, ,0 4 111826 M^ 337931 MjV 



3 



+ 9529520 ^^^^'^ " ~81 M| 1890"!^^ ^ 

6788489 Mi«a;2 _ 13638059 Mfa;^ 
^ 79380 M| 15876 M|~ ^ 

768721 M/o^ 98062 Mfcu^ 20558 Mfu; 
^ 79380 ^^67 M^ ~ 2835 Mq^ ^ 



33 



107608 Mfuj^ 855242 Mfuj'^ 79368645 Af>^ 

105 Mf ^ 405 M| 32032 Uf ^ 
2431 ,10 10 22003 Mf 6435 , o 
256 ° 19845 256 ° 

1155^ ,0,. 4 3003 ^ ,0,. 6 , 4107659 Mf 



128 ° 128 ° 11351340 Mq^ 

218931 M| 131451 323 Mfa;5M| 15999047 Mfcu^Mf 

16016 ^ 5040 Ml 336 M| ^ 

39279559 Mf a;^M| 25118767 Mfu^Ml 
^ 840 Mf 5040 M| ^ 

184244621 Mfuj^M^ 1343113868 Mfu^M'^ 
7560 Ml 315315 M| 

6356431271 Mfu^M^ 21620803 Mjcu^Ms 
183456 1260 Ml ^ 

83616411 MlMiuj^ 25940077 M^Miuj 
^ 28028 Mo 84084 Mq ^ 

q 46612879 Mj^Mscu^ 705291019 M^^Mao; 
+ 662M0M2M1V+ ^ ^ 



2 

+ 



11760 Mo 291060 Mq 

(57) 



^iH = , g2{uj) = -{l-uj^)]^ujMl 

g^{u) = -(1 - u'^) {2>MlMxu'^ - 2Mlu + 3a;MoM2 - Mo^Mi) 

3Mo 



gi{uj) = -{l-uj'^) — 2(126MiVm2 + 147u;2Mo^ + 210Mo^M2u;^ + 

168Afo 

- 882Mia;2MoM2 + 336Mfa;^ - 234Mo^M2u; + lOMo^M^^u; + 

- 63a;Mo^ + I26M1M0M2 - 28Mf ) 

g^{uj) = -{l-Lo'^)^^^{MhMiU}^MlM2 + lSmM^MiUJ^ + h2MlMl 

+ 1260Mo^M2Cj-^ - MQMluj^M^ - 2A<oAM^uj^ - 4032M^M^uj^ + 

+ m^2M2MlMQU^ + 1566Miu;2Mo^M2 - 852M^Mfu^ + 

- 1400Mo^Mitu^ + 1512Mo^M|c<; - 336Mo^M2o; + 588Mi^tu + 

- 2604M2M^MoUJ + 224M^ujM^ + 126MiM^ - 255MiM^M2) 



34 



^eM = -(1 - 55440^4 (^6016Mf + 948640MV - 395472MiV + 

+ 197736Miu^M^M2 - SSieOOMicu^Mg^Ma + lUlUM^MiM^ + 

+ 485056M^uj^M^ - 99616M^uM^ - SSOOSMj^MsMo + 

+ 103950cj^M2^Mo^ + 221760MiVMo^ + 6776M^u;^M^ + 

- 512820MJM2Cj3 - 310410cjM|Mo^ + 785340uj^ + 
+ II319OMJM2CJ + 381150MJM2CJ^ + 533610Mo^MiV + 
+ 78386M^M^u + llOSSMiMg^Ma - U8U0M^uj^M^ + 

- 103950cj^M(}° + 17325cuMo^° + 114345cj^Mo^° - 739271^1^7^0^ + 
+ 488664:M^ujM2M^ - 17 41608 M^uHliMj^ + 

+ 1976436M^M2Uj^Mo - 2029104MiMiUj^M^ + 

- 4162620MfM2u;^Mo + 3735270 M^Miu^M^) 

g,{uj) = -{I- uj^) ^3^^g^^Q^5 (-125789664Mf^^ + 76076000Mf^^ + 

- 7847840Mfa; - 176138820Mf M2u;^Mo^ + 2044328Mi^u;=^Mo^ + 

+ 77882805MjMiu;^M2 - l6l09808Mluj'^ M2MI - 3A&U2M^Ml + 

- 47750703MJMiu;''M2 - l2Alb2mMluM2Ml + 
+ lh26l2mOMluj^ M2MI + 10947807Mia;2MjM2 + 



Al20limM^uj^M^ + 1770912Mi^tuMo^ + 6450444Mi^M2M, 



- 10579590M2MiMo^ + 24594570Mo^Mj^cu^ - 1563962471^0" M2CJ^ 

- 87681880MfcjUfo^ + 32785896Affcj2Mo^ + 179026848M|u;^M3 + 

- 141405264M|cj^Mo3 + 20756736M|cjMo^ - 1017016Aff Af^ + 

- 45939894A^o^Affcj^ + 4926636Afo^cj^2 + 22702680A//o^Af2C^^ + 

+ 45090045Afo^°Aricj^ - 54339285A^o^°A^icj^ + 15786771Afo^°A^icj2 + 

- 15135120Afo^A^i2cj^ - 102162060Afo^cj^A^2^ - 693381AfiAfo^Af2 + 
+ 36159552Afo^o;^Af| - 1153152Afi2tuAfo^ + 12803830Affo;^Afo^ + 

+ 1729728 Afo^Af2C^ + 10426416 Afi^cu^Afo^ - 373258 170Af2^A^ia;^Ai'o^ + 

+ 47411364A4';'Ai'2C^Ai'o + 633008376A^|A^i^^^o + 

- 782581 80 Af2^Afi^o;Afo2 - 435 146712 Af2AfiCU^Afo + 
+ 395417880Afj^Af2tc'^Mo3 + 661368708Af2Mi^cj^Afo + 

- 879458580Af|Af2a;5M2 + 1971555 12Af|Afia;2A^o^ - 675675A^o^°Afi) 
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+ 170270 lOOMfcj^Mo^ - 2115133020Moi^M2CU^ + 
+ 398670272Mfcj^Mo^ - Wdll22m MfuMl^ + 

r8/\#2, ,3 1 '/noi on/io /\/f8 /\#3, ,2 , on'7i fi'/co /\/f8 /\#2, 



- 1343785212M°M|cj^ - 170218048M°Mfcu^ + 89716788M°M|cj + 
+ 14477068Mi^cuMo^ + 175258512M|MiMo^ - 341059716MiVMo6 + 

- 101474100Mo^^M2Cj + 1763421660Mo^M2^cj^ - 12684672MiM^M2 + 
+ 1047642596Mi^cj^Mo6 + 2705402700Mo^°Mi2cj^ + 

- 3787563780Mo^°Mi^cj^ + 121396275Mo^V + 

+ 1229728500Mo"M2Cj^ + 964863900Mo"M2CJ^ + 14240512Mf Mq^ + 

+ 833152320M3^u;2 + 5870184320M3^u;6 + 1040539500Mo^cu^M2 

- 4681316640Mja;^ - 18610592Mi^ - 315630315Mo^^a;^ + 
+ 845644871^0^71^1^ + 225450225M(}^a;^ - 11036025Mo^^a; + 
+ 1929727800MJm2cj^M2 + 3321323688M2M2Mja;3 + 

- 189787224M2Afi2Mja; - 7883403528M2M2mJ6^^ + 

- 1271350080Mia;^Mo^M2 + 103975872Mia;^Mo^M2 + 
+ 400912512Miu;^Mo9M2 - lS2im66AMl Milo'^M^ + 

+ 10867016160M2^Mia;^Mo^ - 2186245776M|Mia;^Mo^ + 

+ 261449760Mfa;M2 - 2945401888Mfu;^Mo2 + 130306176Mf M2M0 + 

- 281465184Mf M|Mo^ + 63608688Mfu;^Mo^ + 5295205344Mfa;^Mo^ + 
+ 181261080M|MiMo^ + 2260177920Mi^u;^Mo^ - 316277104Mfu;^Mo^ + 

- 96782400Mi^M2Mo5 + 6885419688M|u;^Mo5 - 1009570716M|u;Mo^ + 

- 5665579920Mf M2a;^Mo + 30426396000M2Mi^a;^Mo + 
+ 576521088Mf M2a;^Mo^ + 63314747496M|a;^Mf Mq^ + 

- 35939855952M2Mfa;^Mo - 1751809488Mi^M2u;Mo^ + 

- 8770104252M|tu^Mo^ + 27390242880M|a;^MiMo^ + 

+ 11488060584M2^Mfcj^Mo^ + 3265593804M2^M2a;Mo^ + 

- 28341268956M|cj6MiMo3 - 6103020924M|Micu2Mo3 + 

- 57162465360M2^Mfo;^Mo^ - 28659080736M2Mi^a;^Mo3 + 
+ 17840332224M2Mi^cu^Mo3 + 1916785728MfM2Cu2Mo^ + 

- 27760378680M2m2u;^Mo^ + 3975471 1980M|M>5Mo^ + 
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- 10695484800MfM2a;^Mo5) 

+ 5450156712Micu^Mo^4 - 1821910868MiVMo^° + 

+ 2542700160Mo^^M2CU^ + 35239099896M|cj^MJ 

+ 12543230700Mo^°Mfcj^ + 1347869952M|Mo^°c;3 + 

+ 98163208MJMo^ - 4116752640Mo®Mi^cj^ + 

- 876431556Mio;2Mo^^ + 105656346Mj'M2Mj + 

- 327742272M2^Mo^°o; - 55351296M2M(}^o; + 

- 38575496960Mfo;^ - 2924105184Mo^3M2Cc;^ + 
+ 574766192Mi^cj + 5449384512Mo^°M2^cj^ + 

+ 964240992Mo^Mf cj^ - 22746463740Mo^°Mj'cj^ + 

+ 874377504M2Afo^''^cj3 + 38297571312Afi^cj-'^ + 

- 4484584104cuM2MfMo + 81237765MiM^^M2 + 

- 1695133440M(}2^2^^ - 10575665100Mo^^Mia;^ + 
+ 2237664Mo^Mi^o^ + 3009066840M|cuMJ + 

+ 6271615350Mo^^Mia;^ - 16666A3A6 M^M^Ml + 

- 11987015040Mj°M|a;^ + 11804986908Mo^°Mfw^ + 
+ 71001216Mi^a;3Mo^ - 7268072760M|cc;^MJ + 

- 13562878680M|u;^Mj + 19AM03A56M^'^M^u;^ + 

- 9957948000Mfu;^ - 582918336Mi2a;^Mo^^ + 

+ 95484445056M2u;^Mi^Mo^ + 22072050MiMo^^ + 

+ 3621679776Mo9M2M>^ + 10897286A00M^M2M^u;'^ + 

+ 20985789825M(}^Miu;^M2 - 31206455280Mo^^Miu;^M2 + 

+ 48082385880Mo^MiM|a;^ - 79406203500Mo^MiM|a;^ + 

+ 7236479250Mo^MiM2^a;^ - 29285259180Mi^M2a;^Mj + 

+ 992103912Mf M2a;^Mj + 2605132530Mf M2u;^MJ + 

+ 58534088184Mj^M2CU^Mj - 38348741 16Mo^MiM|a;^ + 

- 4918591392M2o^3Mo9M2 - 2773705080MiM(}^M2CU^ + 
+ 532975872Mi2a;Mo^M2 + 15794246754MiAfo"M2u;^ + 

- 772894980Mi^M2Mo^ - 47263087872u;^M2^Mo^ + 



37 



- 5233976176Mlu^M^ - 5796590S0 M^ujM^ + 

- 13021279128Mi^a;^Mo^ - 2494612120Mfa;^Mo^ + 

- 315702232M^u;^M^ - A3A200SS712MJuj'^ + 

- 1443820371M|MiMo^ + 1913406531MfM|Mo^ + 

- 17048199168Mfcj^Afo^ + 5499844064Mfcj''^Mo^ + 
+ 32426263584M/cj^M2 + 5425958984Mfcj^Mo^ + 
+ 1506160656cuM|Mo^ - 18953494560cu^M|Mo^ + 
+ 57315834576cj^M2^Mo^ + 40359368Af(J°Aff + 

+ 4325764872Mi^Af2Cc'Mo^ - 221128005384:M2M^uj^M^ + 

- 226503770769M|cu^MiMo^ + 49793449959M|Micj2Mo^ + 

- 9063952560M2^Mi2cjMo^ + 273714152712M2Mi^cj^Mo^ + 

+ 237115503765M2VMiMo^ + 43921 7769789M2Mj'cu^Mo^ + 

- 500779192221M2VMj^Mo^ - 84523620003M|Mfcu^Mo^ + 
+ 11427580164cuM|MfMo2 + 75509594160M2Mfcu^Mo + 

- 10046350314M2cjM|Mo^ + 269384667552Mfo;^M2Mo + 
+ 38271033372Mf M2Cj2M^ + 6U06m7556M^uj^M^M^ + 
+ 57847286016M^ M^uj^Ml^ + U6m916572M^u;^M^M^ + 

- U0620399920Mluj^ M^M^ + 7UU85480M2MfujHl^ + 

- 480543897834:M^u^ Ml + 418419689688Mi^cj^A'/|M^ + 

- 590669679600 A^i^cj^Af^Afo " 1 83306795672 A^fcj^A^2^Afo + 
+ 4067972064Ar|Afi^^^o - 279791936424A^fa;^A^2Mo + 

- 36119871216Ai'2A4"iCU^Ai'o^) (58) 
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